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Summary. — I briefly review a systematic approximation scheme of QCD in 
which the quark model picture of hadrons emerges at lowest order. A linear A" 
potential arises if Gauss' law is solved with a non-vanishing boundary condition 
at spatial inflnity. Similarly to the Dirac case one can describe relativistic states 
including any number of particle pairs (sea quarks) using valence wave functions, 
whose norms give inclusive probability densities. Provided as(Q^) freezes in the 
infrared, perturbative corrections to the S-matrix can be calculated in the usual 
way, but with states bound by the linear O (Qs") potential instead of plane waves 
in the in and out states. 
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1. Questions and Answers 

Hadrons have two seemingly incompatible features: 



• Their (light) quark constituents are highly relativistic, consequently hadron wave 
functions have abundant contributions from sea quarks and gluons. This is seen 
in the parton distributions measured in hard processes, and is required by the 
underlying theory of Quantum Chromodynamics (QCD). 

• The hadron spectrum reflects just the valence quark {qq, qqq) degrees of freedom. 
Sea quarks and gluons do not manifest themselves as additional states (multi-quark 
states, hybrids, glueballs). The spectra of heavy quarkonia are remarkably similar 
to the atomic spectra of perturbative QED. 



This raises the following questions: 



(*) Presented at the 3rd Workshop on the QCD Structure of the Nucleon (QCD-N'12), 22-26 
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1. What approximation allows relativistic field theory bound states, with prominent 
multi-particle Fock components, to be described by valence quark wave functions? 
The bound states of electrons in a strong Coulomb field have Fock components 
with multiple e^e~ pairs, yet are described by Dirac wave functions that have only 
a single electron degree of freedom. I discuss why this implies that the norm of the 
wave function is an inclusive probability density. I also recall a remarkable feature 
of the Dirac wave functions which has been largely forgotten since the early 1930's. 

2. How does the quark model potential V{r) = or — CpOLs/r arise in QCD? 

The linear term stems from a non- vanishing boundary condition in solving Gauss' 
law for A". The single gluon exchange term is an O (a^) correction to the lin- 
ear term. The perturbative expansion is meaningful provided the strong coupling 
freezes at low momentum transfers, with several estimates indicating as(0) ~ 0.5. 

2. Relevant aspects of the Dirac wave function 

The Born level Green function Gb of an electron in an external A'^{x) field generated 
by a static charge eZ is obtained by summing all Goulomb exchanges between the electron 
and the charge. The residue Rb of a bound state pole, Gb{p^,p) = Rb{E ,p) / [p^ — E), 
satisfies the Dirac equation with potential e^°(a;). In this Born approximation all O (e^) 
vertex corrections as well as loops created by the external field are ignored. Since the field 
is static it does not change the p'^ component of the electron momentum. If p° > —m the 
ie regularization is irrelevant at the negative energy poles of the electron propagators, 
since p° + ^/p^ + rn^ > 0. 

In order to find the wave function of a bound state at an instant of time t we need to 
time-order the interaction vertices. Electrons of negative energy move backward in time, 
creating so called 'Z-diagrams' where one (or several) additional e+e^ pairs are prop- 
agating along with the electron. An equal-time, relativistic bound state wave function 
necessarily contains any number of such pairs. As the strength of the potential increases 
the bound state becomes more relativistic and the pairs become more prominent. 

The usual Dirac wave function ip(t,x) = cxp{—iEt)ip{x) depends on the position 
X oi a single electron, which may have either positive or negative energy. This wave 
function is obtained when retarded, rather than Fcynman electron propagators are used. 
As we noted above, the Green function Gb{p'^,p) is independent of the ie prescription 
at the negative energy poles of the electron propagators. In particular, the bound state 
energies E are the same whether retarded or Feynman electron propagators are used. 
However, the equal-time wave function is obtained after Fourier transforming — > t, 
and is sensitive to the is prescription. In retarded propagation also negative energy 
electrons move forward in time, there are no Z-diagrams and the amplitude of the single 
electron present at time t is given by 'ip{t^ x). 

The fact that the Dirac wave function tp{t, x) is obtained using retarded (rather than 
Feynman) boundary conditions means, in analogy to cross sections [2], that \'ijj{t,x)\'^ 
should be interpreted as an inclusive probability density. As shown by Weinberg [3] , the 
norm of the Dirac wave function should be normalized to unity provided the normalizing 
integral converges. Already in the early 1930's it was realized [4] that the normalization 
integral diverges for most potentials - the 1 /r potential of I? = 3 -I- 1 dimensions being 
an exception. E.g., for the linear potential V{x) = \e'^Z\x\ of QED in D = 1 -I- 1 the 
norm |'0(i,a;)p, and hence the inclusive particle density, is constant at large distances x. 
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This is consistent with the idea of string breaking, particle pairs being created over an 
interval of x where the potential increases by 2m. It is also in accord with the behavior 
of the Dirac wave function for nearly non-relativistic dynamics, e^Z <C m?' [1]. The 
normalization integral being divergent furthermore implies that the Dirac spectrum is 
continuous rather than discrete [4]. Curiously, these properties of the Dirac solutions 
have been scarcely discussed since they were first noticed long ago. 

3. The linear potential 

The linear confining potential of the quark model is necessary for describing the 
hadron spectrum. This introduces a dimensionful constant ~ ^qcd, which is not present 
in the QCD Lagrangian. Since the Schrodinger equation is a Born level approximation the 
constant cannot originate from the renormalization of loop integrals. The only possibility 
consistent with the underlying theory is then to introduce it via a boundary condition. 

In gauge theories the Coulomb potential A'^ is instantaneous, since its time derivative 
does not appear in the Lagrangian. At each instant of time the charge distribution 
determines A'^ through the field equations of motion. E.g., in the Hydrogen atom of 
QED, for each position a;i of the electron and X2 of the proton we find A'^{x;xi,X2) by 
solving Gauss' law, — V^^° = e[S^{x — Xi) — 6^{x — X2)]. The standard solution gives 
rise to the Coulomb potential V{xi — X2) — —a/\xi — X2\ when both the field energy 
^ / d^x(V^°)^ and the interaction energies eA^(xi) — eA^{x2) are taken into account 
(the infinite self-energies are independent of Xi,X2 and thus irrelevant). 

The —a/r Coulomb potential of the Hydrogen atom results when Gauss' law is solved 
with the boundary condition that A'^{x;xi,X2) — > as \x\ — >■ 00. If we instead require 
that the field strength approaches a non-vanishing constant, (VA")^ as |a;| c>o, 

then we have to add a homogeneous solution to Gauss' law, A\ = £-x, where £{xi, X2) 
is a unit vector. The square of contributes a term oc K'^V to the field energy {V 

is the volume of space). This term is irrelevant provided A is a universal constant, 
independent of Xi,X2. The field energy arising from the interference of A\ with the 

(e) potential is finite provided the state is neutral, and is cx eK t ■ {xi — X2). Choosing 

1 \\ Xi — X2 makes the action stationary under variations of the unit vector £, preserves 
rotational invariance and gives rise to the linear potential Va oc eA^|a;i — X2\. 

A linear potential arises also in QCD through an analogous boundary condition [5]. 
This does not explain why we should choose A = in QED and A ~ Aqcd 7^ in 
QCD to describe observed phenomena. However, it shows that a linear potential is 
compatible with the field equations of motion. A is a hidden parameter, not present in 
the Lagrangian, whose value must be determined by experiment. As indicated above, 
the solution requires A to be a universal constant and the state to be charge neutral. In 
QCD bound state solutions were found [5] for color neutral qq mesons and qqq baryons. 
No analogous solutions exist for states with a higher number of valence quarks. 

4. Born level bound states 

The linear A^ potential ensures confinement and relativistic binding. Since the po- 
tential is instantaneous hadrons may be described using "inclusive" valence quark wave 
functions analogous to those of Dirac states. Perturbative corrections can then be added 
systematically. In the S'-matrix expression Sfi ~ out(/|T[cxp(— i J d^xT-Lj)^ the in 
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and out states should consist of neutral states bound by the linear potential, instead of 
the free plane wave states of standard perturbation theory. In effect, one expands around 
lowest order states which already have the basic features of hadrons such as color con- 
finement. As in the Taylor expansion of ordinary functions, the full perturbative series 
formally gives the exact result independently of the starting configuration. 

For the perturbative expansion to be meaningful the strong coupling should freeze at 
a perturbative value in the infrared. Phenomenological and theoretical studies [6] find 
as(0) ~ 0.5, which justifies the use of perturbation theory. 

The Poincare invariance of the S'-niatrix requires that the in and out states transform 
covariantly. Equal-time, interacting states transform dynamically under boosts, since the 
concept of equal time is frame dependent. The Poincare invariance of the action, and 
the Lorentz invariant formulation F^^F^^'^ = —2 A'' of the boundary condition, suggests 
that the bound states should transform correctly under boosts. A detailed analysis in 
D = 1 -I- 1 dimensions showed that the boost operator indeed transforms the states as 
required. From this follows, in particular, that the bound state energy E has the correct 
dependence on the CM momentum, E = V + P^. The same dependence was found 
previously in Z? = 3 -I- 1 [7]. The covariance holds only for an exactly linear potential. 
Under boosts the D = 1 + 1 wave fimction contracts oc 1/[E — V{x)]. In the non- 
relativistic limit y <C -B this reduces to the usual Lorentz 7-factor, but in general the 
wave function contracts at a rate which depends on the separation x of the constituents. 
In effect, the relevant quantity is thus found to be the kinetic energy p° — rather 
than the canonical energy p'^ familiar from classical physics. 

For a more complete list of references and a discussion of the wave functions, form 
factors and parton distributions of the bound states in Z) = 1 + 1 see [1] . 
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